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Spin-orbital entanglement in the ground state of a one-dimensional SU(2)®SU(2) spin-orbital model is 
analyzed using exact diagonalization of finite chains. For S — 1/2 spins and T = 1/2 pseudospins one 
finds that the quantum entanglement is similar at the SU(4) symmetry point and in the spin-orbital valence 
bond state. We also show that quantum transitions in spin-orbital models turn out to be continuous under 
certain circumstances, in constrast to the discontinuous transitions in spin models with SU(2) symmetry. 
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Rich magnetic phase diagrams of transition metal oxides and the existence of quite complex magnetic 
order with coexisting ferromagnetic (FM) and antiferromagnetic (AF) interactions, such as A-AF phase in 
LaMnO,3 or C-AF phase in LaVC>3, originate from the intricate interplay between spin and orbital degrees 
of freedom — alternating orbital (AO) order supports FM interactions, whereas ferro orbital (FO) order 
supports AF ones [1|. While in many cases the spin and orbital dynamics are independent from each 
other and such classical concepts apply, the quantum fluctuations are a priori enhanced due to a potential 
possibility of joint spin-orbital fluctuations, particularly in the vicinity of quantum phase transitions J2j- 
Such fluctuations are even much stronger in t2 g than in e g systems and may dominate the magnetic and 
orbital correlations [3], which could then contradict the above classical expectations in certain regimes of 
parameters. Recently it has been realized |4| that this novel quantum behavior is accompanied by spin- 
orbital entanglement, similar to that being currently under investigation in spin models Q. 

In general, any spin-orbital superexchange model derived for a transition metal compound with a per- 
ovskite lattice may be written in the following form: 



n 

7 



^orb, (1) 



where 7 = a,b,c labels the cubic axes — depending on the direction of a bond (ij) the interactions 
take a different form. The first term in Eq. (Q~|) describes the superexchange interactions ( J — 4t 2 /U 
is the superexchange constant, where t is the hopping element and U stands for the Coulomb element) 

"(fy) " (~y) 

between transition metal ions in the d n configuration with spin S. The orbital operators J^ ' and K^J 
depend on Hund's exchange parameter rj = Jh/U, which determines the excitation spectra after a virtual 
d™d™ —> d™ +1 d™~ 1 charge excitation. Therefore, realistic models of this type (some examples were given 
recently in Ref. [6] and analyzed using the mean field approximation in Ref. Q) are rather involved and 
contain several terms). In addition, orbital interactions can also be induced by the coupling to the lattice, 
and appear in H ib term which depends on a second parameter V. 

In the limit of r\ = 0, however, and for t2 g orbitals, and operators simplify and contain only 
a scalar product Tj • Tj of T = 1/2 pseudospin operators which stand for two active t^g orbitals along 
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Fig. 1 Two degenerate SOVB ground states |$±) 
with alternating spin (S = 0) and orbital (T 
(I) singlets for the SU(2)xSU(2) model 101 at the 
dimer point p = 3/4. 



7, so the SU(2) symmetry follows both for spin and for pseudospin interactions. Here we shall discuss 
primarily a one-dimensional (ID) model in this idealized situation and investigate the dependence of spin- 
orbital entanglement on the type of underlying interactions. To characterize the ground state we evaluated 
intersite spin, orbital and composite spin-orbital correlations, defined as follows for a bond (ij) [4|: 



Sy = 0i ■ S j )/{2Sf, 

Cij = [(0i ■ • 7>)) - (Si ■ ■ %)] /(2S) 
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Note that quantifies the quantum entanglement — if Cy < spin and orbital operators are entangled 
and mean field approximation, i.e., decoupling of spin and pseudospin operators is not justified. 
The ID SU(2)xSU(2) model, 



Hj = jJ2 (3 ' + x ) ' + y) > 



(4) 



has two parameters x and y. Its phase diagram in the (x, y) plane consists of five distinct phases which 
result from the competition between effective AF and FM spin (AO and FO pseudospin) interactions on 
the bonds [ 8 1 . First of all, the spin and pseudospin correlations are FM-FO, Sij = Tij = \, if x < — \ and 
y < — j. Then the ground state is characterized by the maximal values of both total quantum numbers, 
S = T = N/2, where N is the chain length, its degeneracy is d = (N + l) 2 , and the quantum fluctuations 
are suppressed. Two other phases, with either spin FM or pseudospin FO correlations, show also no 
entanglement as the wave function factorizes, i.e., spin and orbital operators may be then still decoupled 
from each other. In these cases only either spin or pseudospin quantum fluctuations occur. On the contrary, 
in the other two phases both spin and orbital correlations are negative [9|, so one expects that quantum 
entanglement takes over. 

This motivates our study of the SU(2)<g)SU(2) model along the line p = x = y, where S — 1/2 spins 
and T — 1/2 pseudospins appear on equal footing. However, in order to allow also for a larger value of 
spin S per site, we write the Hamiltonian along this line in the following way, 



J 



S l+1 + - P S{S+l) )(Ti-T i+1 + P 



(5) 



LU 




Fig. 2 Energy spectrum for the S = 
1 /2 spin-orbital chain {5} of TV = 4 
sites for increasing p. Some of the 
256 eigenenergies have high degen- 
eracy d. At p = —0.25 the ground 
state changes from high spin-orbital 
state (S = T — 2, d = 25) to spin- 
orbital singlet state (S — T = 0, 
d — 1); its degeneracy atp = —0.25 
is d = 194. The SU(4) point (p = 
0.25, d = 1) and the SOVB point 
(p = 0.75, d — 2) are marked by 
vertical dashed lines. 
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Fig. 3 Intersite spin and orbital correlations Sij = Tij © (circles), and composite spin-orbital correlations dj ® 
(crosses), as obtained for the ID spin-orbital chains for increasing p in Eq. (5J> with: (a) N = 4, and (b) N = 8 sites. 



with ^pS(S + 1) = p for S = 1/2. Interestingly, for p = | the model given by Eq. © has an exact 
doubly degenerate ground state |$±), with alternating spin and orbital singlets forming a spin-orbital 
valence bond (SOVB) phase (Fig. 1). These states resemble the Majumdar-Ghosh valence bond states 
in the ID spin model with next-nearest interactions [10] and have exact energy E = given by either 
spin or orbital singlet. Each |$±) state is a matrix product state in both spin and orbital sector JTTJ . For 
5=1/2 and p = \ one recovers the celebrated SU(4) model, with all three correlation functions: £y, 
T^, and • Si+i)(Tj • Tj + i)) being equal to each other fl2l . The energy spectra in these two cases 

look quite differently for an N — 4 site chain — an equidistant energy spectrum with high degeneracies of 
each excited state with energies (—1 + J where n = 0, 1, • • • , 4 is found at the SU(4) point, while the 
spectrum consists of many eigenenergies with lower degeneracies at the SOVB point (Fig. 2). It is a priori 
not clear in which of these two points the quantum entanglement is stronger, but one might expect that it 
would be more pronounced at the SU(4) point. 

In fact, it was shown recently that a reduced von Neumann entropy is there maximal [13]. While the 
spin-orbital entropy could help to understand the consequences of the entanglement at finite temperature 
which might trigger phase transistions to phases with strong dimer correlations [ 14] observed in experiment 
1 15 1, here we suggest that a useful measure of entanglement in the ground state is the intersite spin-orbital 
correlation function © which quantifies the error of the mean field decoupling of spin and orbital operators 
on individual bonds. The intersite spin, orbital and spin-orbital correlations demonstrate a quantum phase 
transition between the high spin-orbital FM-FO state (S = T = 2) and the singlet entangled state (S = 
T = 0) at p = — j in the ID model © for N = 4 and N = 8 sites, see Fig. 3. In fact, the Hamiltonian 
(© is then a product of singlet projection operators in spin and orbital space, so the FM-FO state has the 
lowest possible energy Eo = and is degenerate with several other states (Fig. 2). 

As expected, dj — for p < — \, and the spin and orbital degrees of freedom disentangle. The 
situation is qualitatively different above the quantum transition at p — j, where Cy ~ —0.25 (CV, = 
—0.2595, —0.25 for an TV = 8, 4 site chain at p = —0.249, respectively). Although we are not interested 
here in accurate quantitative values of the respective intersite correlations, we observe that the N = 8 
site chain comes already quite close to the thermodynamic limit at the SU(4) point — we found Sij — 
— —0.223 in place of —0.215 for a chain of N = 100 sites (which almost reproduces the Bethe ansatz 
result for an infinite chain) fl2l . Furthermore, in the case of N = 8 site chain, the data show a smooth 
decrease of Sij = TV, correlations with increasing p, indicating the tendency towards AF-AO spin-orbital 
correlations in the regime of large p. This is markedly different from the N = 4 site chain, where these 
correlations are constant (Sij = T^ = —0.25) in the entire range of — | < p < |, i.e., up to the SOVB 
point, where the data suggest a second quantum phase transition to the AF-AO state (Fig. 3a). In fact, this 
transition is a finite size effect and is replaced by a smooth crossover towards smaller values of both Sij 
and Tij correlations in an N = 8 site chain (Fig. 3b). 
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Fig. 4 Ground state energy Eo (diamonds), and mean field energy Emf (squares) per bond, as obtained for the ID 
spin-orbital chains for increasing p in Eq. l[5}, with: (a) N = 4, and (b) N — 8 sites. 

It is quite remarkable that the composite spin-orbital correlation function (0l has a shallow minimum 
Cij — —0.2812 for an N = 8 site chain at the SOVB point, in spite of decreased individual spin and orbital 
correlations = Ty = —0.3749, suggesting rather independent spin and orbital dynamics. Therefore, 
precisely at this point the entanglement in the exact wave functions |$±) shown in Fig. 1 is larger than at 
the SU(4) symmetric point, where CV, = —0.2667. We emphasize that sufficiently long chain with N = 8 
sites was necessary to reach this conclusion, while shorther chains give inconclusive results, either due to 
particular stability of the SU(4) singlet for N — 4 sites, or due to frustrated four-site singlet correlations 
for N — 6 sites which result in reduced values of |C,-j | and Emf (Table I). 

In the entire regime of singlet states (S = T = 0) in the ID spin-orbital model Eq. (O, one finds large 
corrections to the mean field energy normalized per one bond, 

Emf = ((& • +p) ((3 ■ fi+i) +p), (6) 

as shown in Fig. 4. This demonstrates that one should not decouple spin and orbital operators from each 
other as this procedure suppresses an essential part of joint spin-orbital quantum fluctuations and leads to 
inconclusive results. For the same reason, effective magnetic exchange constants obtained by averaging 
over orbital operators do not represent a useful concept and cannot be introduced in the entangled regime, 
similar to the realistic spin-orbital models for transition metal oxides J6l|7). 

We have verified that quantum entanglement is significant in the SOVB wave functions (Fig. 1) also 
at higher values of spin S in Eq. (0. However, when spins increase, the AF correlations weaken and AO 
correlations are enhanced. This may be seen as partial decoupling of slow and fast quantum fluctuations 
associated with large (small) spin (pseudospin) value, respectively. Indeed, the values of \Cij | and of -Emf 
decrease with increasing S, as shown in Table I. 

Table 1 Intersite spin correlations Sij, orbital correlations Tij 10, and composite spin-orbital correlations 
dj ®, as well as mean field -Emf energy per site (|6}, obtained for various SOVB phases: for S = 1/2 
with increasing cluster size of N = 4, 6 and 8 sites, and for S = 1 and S = 3/2 with clusters of N — 4 
sites. By construction, the exact ground state energy is Eo — in all cases. 



s 


N 


Sij 


Tij 




Emf ( J) 


1/2 


4 


-0.3500 


-0.3500 


-0.1600 


0.1600 




6 


-0.3735 


-0.3735 


-0.1417 


0.1417 




8 


-0.3749 


-0.3749 


-0.2812 


0.2812 


1 


4 


-0.2429 


-0.3643 


-0.0992 


0.0992 


3/2 


4 


-0.2050 


-0.3690 


-0.0806 


0.0806 
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Fig. 5 Intersite spin Sij (filled circles), orbital Ty (empty circles), and composite spin-orbital dj (crosses) corre- 
lations, as obtained with an N = 4 site chain for increasing Hund's exchange r\ in: (a) the vanadate model (TJ of 
Ref. [H, and for (b) the same model without orbital fluctuating terms in Eq. {T}. Orbital interactions resulting from 
GdFeC>3-type distortions are given by 7i or b = — V J^. T*Tl +1 with V = J, and favor FO order at small rj. 



Finally, we would like to emphasize that the quantum entanglement is also a common feature of spin- 
orbital models (1) derived for transition metal oxides when the parameter rj = Jh/U is small |4). (Note 
that large rj plays here a similar role to small p in the SU(2)®SU(2) spin-orbital model.) Also in such 
cases mean field procedure fails and joint spin-orbital fluctuations Cij dominate in the ground state. Such 
models exhibit even richer behavior than the idealized SU(2)®SU(2) spin-orbital model considered above. 
In fact, the SU(2) symmetry concerns only spin interactions while it is removed in the orbital sector due 
to the analytic structure of Coulomb interactions and the multiplet structure of excited states at 77 > 0. As 
a result, one may find continuous orbital transitions of the crossover type even when the orbital quantum 
number T changes, as shown in Fig. 5a for the vanadate model of Ref. 1 14| with 5=1 spins. However, 
when the respective interactions are simplified to classical Ising terms (Fig. 5b), such transitions appear to 
be first order JT6), and are qualitatively similar to those encountered usually in the spin sector. 

Summarizing, we have demonstrated large quantum entanglement in the SU(2)®SU(2) spin-orbital 
model for S = T = 1/2, which appears to be somewhat stronger in the exact spin-orbital valence bond 
states of Fig. 1 than at the SU(4) symmetric point. While the spin-orbital entanglement occurs in a 
discontinuous way in quantum spin transitions, it may appear gradually in continuous quantum transitions 
which involve orbital degrees of freedom. 
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